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FEEDBACK  CONTROL  FOR  AN  ABSTRACT 
PARABOLIC  EQUATION 

Rouben  Rostamian1  and  Thomas  I.  Seidman2 

Department  of  Mathematics 
University  of  Maryland  Baltimore  County 
Catonsville,  Maryland  21228,  U.S.A. 

Takao  Nambu3 

Department  of  Mathematics 
Kumamoto  University 
Kumamoto,  Japan 


1.  Introduction 

In  our  previous  paper  [6]  we  considered  the  problem  of  the  exponential  stabiliza¬ 
tion  of  the  heat  equation  with  Neumann  boundary  conditions 

1This  research  has  been  partially  supported  by  the  Office  of  Naval  Research  and  by  the  National 
Science  Foundation  under  grants  N00014-86-K-0498  and  DMS-86-02006  respectively. 

2  This  research  has  been  partially  supported  by  the  Air  Force  Office  of  Scientific  Research  and  by  the 
National  Science  Foundation  under  grants  AFOSR-82-0271  and  CDR  85-00108,  respectively.  Part 
of  this  work  was  done  while  this  author  was  visiting  at  the  Systems  Research  Center,  University  of 
Maryland. 

3  This  research  has  been  partially  supported  by  the  Air  Force  Office  of  Scientific  Research  under 
grant  AFOSR-82-0271  while  this  author  was  visiting  the  Department  of  Mathematics,  UMBC. 


(1.1) 


in  n  x  (0,  oo) 
on  <9fl  x  (0,  oo) 
x  in  H 


2 


du 

dt 

du 


—  A  u 


„  =  0  on  dQ  x  (0.  oo) 

du 

u(x, 0)  =  u0(x)  x  in  H 

in  the  smooth  and  bounded  domain  fi  in  Rn,  where  d/du  denotes  the  differentiation 
in  the  direction  normal  to  the  boundary  <90.  The  solutions  of  this  problem,  although 
stable  (in  the  sense  of  Liapunov,  considered  for  any  Lp  space,  say)  are  not  asymp¬ 
totically  stable,  because  spatially  homogeneous  functions  do  not  decay.  A  practical 
way  to  induce  asymptotic  stability  is  to  introduce  a  feedback  mechanism  which  ob¬ 
serves  the  temperature  at  a  patch  of  the  boundary  dO  and  causes  an  appropriate 
exchange  of  heat  through  the  boundary  or  through  an  interior  subset  of  O.  A  feed¬ 
back  mechanism  of  the  latter  type  may  be  modeled  by  introducing  a  term  of  the 
form 


(1.2)  $(M)  ~  -£  (fgnu(y>t)<P(y)dy^  er(x) 

on  the  right  hand  side  of  the  equation  in  (1.1).  The  functions  ip  and  a ,  prescribed 
on  <90  and  O,  respectively,  are  indicators  of  the  sites  where  the  observation  and  the 
feedback  take  place.  The  coefficient  e  is  the  ‘gain  factor’  of  the  feedback.  Feedback 
mechanisms  of  the  former  type  may  also  be  treated  by  combining  the  method  of 
section  6  of  [6]  and  the  estimates  developed  in  the  present  paper. 

More  generally,  we  may  consider  a  more  complex  diffusion  processs  governed  by 
a  general  (autonomous)  linear  elliptic  operator  A,  and  several  feedback  mechanisms 
of  the  type  above  that  act  simultaneously  and  independently,  so  that  the  governing 
equation  becomes 


(1.3) 

For  instance, 


du 

dt 


Au  (Jgn«MMy))  ai(x)- 


n  d 

Au=  V  4- 
u=l  dx‘ 


Otj  (x) 


du 

dxj 


shut  \du 


may  correspomd  to  a  heat  conduction  problem  in  an  anisotropic  and  inhomogeneous 
medium. 


The  method  introduced  in  [6]  for  the  stability  analysis  extends  to  this  and  more 
general  situations.  In  this  paper  we  study  an  abstract  evolution  equation  which 
includes  the  problem  in  [6]  and  the  problem  corresponding  to  equation  (1.3)  as  special 
cases.  Our  main  result,  in  the  context  of  the  latter  problem  is,  roughly,  that  if  certain 
hypotheses  are  satisfied  and  if  e  is  sufficiently  small  then  the  system  is  exponentially 
asymptotically  stable.  Studies  of  similar  kind,  albeit  with  quite  more  complicated 
sufficient  stability  conditions,  appear  in  various  places  in  the  literature;  we  refer  the 
reader  to  the  articles  in  [l],  [5],  [7],  and  [8]. 

The  smallness  of  the  feedback  gain  factor  e  is  essential  for  stability.  In  [6]  we  give 
an  example  where  for  a  large  e  the  solution  grows  exponentially  as  t  — >  oo.  For  a 
comprehensive  study  of  the  spectral  properties  and  the  spectral  resolution  of  elliptic 
operators  of  the  type  occuring  in  (1.3)  see  the  recent  article  of  van  Harten[3j. 

2.  The  Abstract  Problem 

Let  A  be  a  sectorial  operator  (cf.  [4])  on  a  Hilbert  space  H,  and  denote  by  etA 
the  analytic  semigroup  generated  by  A.  Suppose  that  the  null-space  Hq  of  A  is 
finite-dimensional  and  non-trivial,  and  that  the  following  hypotheses  hold: 

Hypotheses  on  A: 

(i)  There  exists  a  continuous  projection  operator  Q,  not  necessarily  orthogonal, 
of  H  onto  Ho,  which  commutes  with  A\  thus  AQ  =  QA  =  0. 

(ii)  If  P  =  I  —  Q,  where  I  is  the  identity  mapping  in  H ,  and  Hi  =  range(P)  = 
null  space(Q),  then  there  exist  positive  constants  a  and  a  such  that 

(2.1)  etAv  <  ae~at\v\,  v  E  Hi. 

Thus,  in  the  dynamical  system  corresponding  to  etA,  all  orbits  in  Hi  are  asymp¬ 
totically  exponentially  stable.  The  points  in  H0,  on  the  other  hand,  are  equilibria  and 
remain  stationary.  We  now  introduce  a  feedback  process,  motivated  by  the  concrete 
application  in  [6],  which  makes  all  orbits  in  H  tend  exponentially  to  zero  as  t  — »  oo. 


We  suppose  we  are  given  q  elements  {ai,  -  •  •  ,aq}  of  H  and  p  linear  (possibly 
unbounded)  functionals  {/1} •  •  •  ,lp}  on  H,  for  which  there  exists  a  p  X  q  real  matrix 
such  that  the  feedback  operator  B,  formally  defined  by 

(2.2)  = 

fc=l j= l 

satisfies  the  following  hypotheses: 

Hypotheses  on  B: 

(i)  For  each  t  >  0  and  k  =  1,2,  •  •  •  ,p,  the  expression  Ik  (e*A)  defines  a  bounded 
linear  functional  on  H. 

(ii)  There  exists  a  <5  E  (0,  l)  such  that 

(2.3)  | lk  (etAu)  |  <  at~s e~at\u\,  u  E  Hu 
and 

(2.4)  \k{u)\  <  a\u\,  u  e  H0. 

(iii)  If  we  define  the  operator  B0  :  H0  —*■  H0  by 

p  i 

(2.5)  B0w  =  -  YlYl"Uklk{w)Qoj,  w  G  H0, 

k= 1  j= 1 

then  there  exist  positive  constants  b  and  /?  such  that 

(2.6)  etB°u;j  <  6e_^t|w|,  w  E  Ho. 

Remark :  The  inequalities  (2.4)  and  (2.6)  imply  the  estimate 

(2.7)  |/*  (etB°ivj  |  <  abe~^\w\,  w  E  Ho. 

Note  that,  formally,  Bo  =  QBQ,  and  that  the  operator  BetA  is  bounded  on  H  for 
each  t  >  0.  In  view  of  this  we  are  able  to  show  that  the  equation 

u  =  Au  +  eBu 

(2.8) 

«(0)  =  Uq 


makes  sense  in  H.  We  prove  the 
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Stabilization  Theorem:  There  exists  a  number  eo  >  0,  depending  on  the  constants 
a,a,b,/3,6  entering  in  estimates  in  (2.1),  (2.3),  and  (2.6)  as  well  as  on  p,q,^jk, 
and  the  norms  of  the  aj’s,  such  that,  for  any  e  G  (0,  e0),  the  solution  of  (2.8)  decays 
exponentially  as  t  —*■  oo. 

We  defer  the  proof  to  the  next  section  but  first  make  some  remarks: 

(i)  It  is  possible  to  obtain  an  explicit  upper  bound  eo  in  terms  of  the  quanti¬ 
ties  mentioned  above.  The  computation  proceeds  along  the  lines  described  in 
section  4  of  [6]. 

(ii)  It  is  worthwhile  to  observe  that  in  our  computation  we  use  the  estimates  in 
(2.1),  (2.3),  (2.4),  (2.6)  very  sparingly.  In  fact,  the  bulk  of  the  computation  in 
the  following  section,  up  to  the  derivation  of  the  system  of  integral  equations 
(3.7),  is  in  terms  of  equalities  rather  than  estimations.  This  has  the  advantage 
of  keeping  the  results  as  sharp  as  possible;  moreover,  the  general  theorem  on 
integral  equations  developed  in  [6]  applies  here  directly. 

(iii)  The  inequality  in  (2.1)  applies  in  situations  where  the  initial  condition  is 
in  the  space  H.  This  assumption  is  of  course  unnecessarily  strong;  in  concrete 
cases  the  initial  conditions  are  often  allowed  to  be  outside  H.  For  instance, 
H  may  be  L2( fi)  whereas  u0  is  a  distribution.  The  regularizing  effect  of  the 
analytic  semigroup  puts  the  orbit  in  H  in  positive  times.  The  estimate  (2.1) 
can  be  modified  to  apply  to  such  irregular  data;  see  [6]  for  details. 

(iv)  The  estimate  in  (2.6)  is  an  assumption  of  the  positivity  of  the  operator  —  Bo 
on  the  space  Ho.  One  might  also  consider  a  possible  order-preserving  property 
of  —B0,  in  the  form 


w  >  0  =>■  —Bqw  >  0,  Vu;  G  Hq 


where  the  inequalities  are  understood  in  the  sense  of  a  partial  order  on  H  for 
which  one  has  a  maximum  principle  for  A.  It  seems  likely,  but  we  have  not 
verified  this,  that  under  this  stronger  condition  the  equation  (2.8)  is  asymptot¬ 
ically  stable  with  no  restrictions  on  the  magnitude  of  the  coefficient  e.  See  [6] 
for  further  discussion  and  practical  impications. 

(v)  It  is  possible  to  show,  although  we  do  not  do  it  here,  that  the  operator 
A  +  eB  entering  in  (2.8)  is  sectorial  for  arbitrary  e.  Hence  it  generates  an 
analytic  semigroup  on  H. 

3.  The  Proof  of  the  Stabilization  Result 

We  apply  the  projections  P  and  Q  to  the  evolution  equation  (2.8)  and  use  the 

notation  v  —  Pu,  w  =  Qu,  to  arrive  at  the  following  coupled  system  of  differential 

equations: 

p  <2 

(3.1)  v  =  Av  -  eJ^JZljkhi^Paj 

k=ij=i 
p  i 

(3.2)  ib  =  -^J2Yl^jklk{u)Qaj 

k=lj=l 

We  also  set  i>o  =  u(0)  =  Pu( 0)  and  wq  =  w(0)  =  Q«(0).  Note  that  the  equation 

(3.2)  may  be  re-written,  in  view  of  the  definition  of  the  operator  B0  in  (2.5)  and  the 

identity  u  =  v  +  w,  as: 

p  ? 

w  +  eB0w  =  ljklk{v)Qoj. 

k= 1 j= 1 

Applying  the  variation  of  constants  formula,  we  get 

(3.3)  v(t )  =  etAv0  ~  l  \e^~a)APaj\  lk{v{s )  +  w(s))  ds, 

k= 1 j=l 

(3.4)  w(t)  =  eetBow0  ~  f  ^jk  [ee(t_s)Bo<5cryl  lk(v(s ))  ds. 

k=  1  j=l 

Now  define 


Pi{t)  =  U{v(t)), 
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A ik(t)  =  ^2  likli  (etAP CTj'j  , 

3= 1 
q 

Vik{t)  =  Y^ljkh  [eetB°Qo2)  , 

i=i 

and  apply  /,•  to  both  of  equations  (3.3)  and  (3.4)  to  obtain 

P  w 

Pi{t)  =  U  (eiAvo)  ~eYl  -M*  -  5)  (p4s)  +  J*M5))) 

jfc=:r0 

li[w(t))  =  /i(e£f%o)-fE  f  »7ii(i--s)pfc(3)ds. 

Jk=i 

Substitute  the  second  of  these  equations  in  the  first,  interchange  the  order  of  inte¬ 
gration  in  the  resulting  double-integral  (justified  by  the  estimates  in  (2.1)  and  (2.3),) 
and  simplify,  to  arrive  at 

Pi(t)  =  U  (etAv2)  JQ  *<*(*  ~  s)lk  ( e£sBw0 )  ds 

k~l 

P  rt 

Xik{t  -  s)pk(s)ds 

k= 1 Jo 

P  rt  rs 

(3.5)  +e  zZ  /  /  ^ik{t-s)rikv{s-T)pk'{T)dTds. 

k,k'=iJo  Jo 

Further  simplification  is  achieved  by  letting  p*  (f )  denote  the  sum  of  the  first  two 
terms  on  the  right  hand  side  of  (3.5)  (i.e.,  the  terms  which  contain  the  initial  data) 
and  setting 

p  rt 

(3.6)  Hij(t)  =  A ij(t)  Kk{r)Vkj{t  ~  r)  dr. 

k= i  Jo 

After  a  bit  of  manipulation,  (3.5)  then  yields  the  following  system  of  integral  equa¬ 
tions: 

(3.7)  Pi(t)  =  p* (t)  -  e  [  »ij{t  -  T)pj{T)  dr 

j= i Jo 

for  i  =  1,2, •  •  •  ,p. 

It  is  easy  to  see  that  by  the  assumption  0  <  6  <  1  in  (2.3),  the  functions  p*  and 
have  integrable  singularities,  hence  the  system  of  Volterra  equations  (3.7)  has  a 
unique  solution.  By  a  computation  similar  to  the  one  in  [6],  it  is  possible  to  show 
that  this  solution  decays  exponentially  as  t  — >  oo,  if  e  >  0  is  sufficiently  small.  Going 


back  to  (3.4),  recalling  that  lk(v(t))  =  pk{t),  and  using  the  estimations  in  (2.1)  and 
(2.3),  we  conclude  that  the  norm  of  w(t)  decays.  Now  since 

lk{u{t ))  =  h(v(t))  +  lk(w(t))  =  pk(t )  +  /*(«>(*)) 

we  get  exponential  decay  for  lk[u[t)).  Finally,  using  this  in  (3.3)  we  obtain  exponential 
decay  for  v(t)  and,  consequently,  for  u(t)  =  v(t)  +  w(t).  This  completes  the  proof  of 
the  stabilization  result. 


Remark :  A  variant  of  the  control  problem  treated  in  this  paper  is  that  of  the  bound¬ 
ary  observation  and  boundary  control.  In  this  situation  the  feedback  mechanism 
involves  interaction  with  fl  entirely  through  the  boundary  of  0.  The  corresponding 
mathematical  model  is: 


ir  =  £  <*(*) 


=  -'Em  o, 


u(x,  0)  =  u0(x) 


in  H  x  (0,  oo) 
on  dfl  X  (0,  oo) 
x  en. 


Abstractly,  this  corresponds  to  having  etAOj  £  H  for  t  >  0  rather  than  Oj  £  H; 

we  must  then  assume  an  estimate  like  (2.3)  for  etAOj  ,  corresponding  to  having  Oj  £ 

D(AS).  It  is  possible  to  show  that  this  system  is  also  asymptotically  exponentially 

stable  provided  that  the  coefficient  e  is  sufficiently  small.  The  proof  is  along  the  lines 

of  the  argument  in  section  6  of  [6],  combined  with  the  generalized  frmework  that  we 

have  developed  in  the  resent  paper. 
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